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The current correlator method has been shown to be a practical tool to extract the charm quark 
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proach using experimental electron-positron annihilation cross section data. We report on the 
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1. Introduction 



The charm quark mass and strong coupling constant are fundamental parameters of the Stan- 
dard Model and thus there is an interest per se in their calculation. They are essential input param- 
eters for the calculation of processes involving charm quarks, such as inclusive radiative B-decays 
and exclusive Kaon-decays [1]. Moreover they play an important role in the estimation of CKM 
matrix elements and the search for new physics beyond the Standard Model [2]. 

Recently the HPQCD collaboration extracted the MS charm quark mass and strong coupling 
constant using temporal moments of charmed lattice current correlators [3]. Using the highly 
improved staggered quark action and a Bayesian prior fitting analysis a few percent precision could 
be reached. Here we report on an ongoing effort to apply this method using a different fermion 
discretization, namely Wilson twisted mass Lattice QCD. In this work, we will not rely on any 
Bayesian prior in the fits of our data and it is one of our goals to understand, whether a similar 
accuracy can be reached as given in [3]. 



2. Low momentum expansion of polarization functions in perturbative QCD 

The general strategy of the current correlator method is the non-perturbative estimation of 
derivatives M n of the polarization functions of in our case the pseudoscalar and vector currents 
from lattice data and to compare them to their continuum counterparts determined in perturbation 
theory. The derivatives are readily deduced from the momentum expansion of the polarization 
functions in the limit q 2 <C m 2 

q 2 U K c = ij d 4 xe iqx (0\T{J K (x)J K (0)}\0), 
{-qhnv+q^v) n c 5 + qnq v n 5 cL = ij d 4 x^ x (0\T{jS (x)4(0)}\0) , (2.1) 

€f = L (^yc^LJm^)), (2.2) 



m>0 



(2.3) 



q 2 =0 



with S = v,a, K = p,s, J p = cjsc, J s = cc, = cy^c, = cy^jsc. The perturbative expansion 
of the coefficients C*' S has nowadays reached the 4 loop level (^ (a v 3 )) ( cf. [4] and references 
therein). 



3. Lattice Formulation 

The calculation we report on here is based on gauge configurations produced by the Euro- 
pean Twisted Mass collaboration (ETMC) using Nf = 2 flavors of maximally twisted and mass 
degenerate Wilson fermions. We refer the reader to ref. [5] and references therein. We treat the 
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charm degrees of freedom in a partially quenched framework by adding a doublet of heavy quarks 
X = (X+ > X ) i n tne valence sector with valence quark action 

^vai = Y,X(x) (Dw + m CI + in h Y 5 T 3 )x(x). (3.1) 

x 

In this framework automatic & (a) improvement [6] is in place with the same critical Wilson mass 
m cr as used in the light quark sector. For the lattice operators representing the physical charm 
currents for a given spin structure F we have three natural choices (given in the physical basis). 

= y/r®li|A; / r = #T<g)T 3 v; Jr = W^^W- ( 3 - 2 ) 

At non-zero lattice spacing the correlation of these operators with themselves will give a different 
results for each operator due to lattice artifacts. Concerning the physical charm fields we would 
need to use the singlet currents and their corresponding translation in terms of the x fields in 
the twisted basis. However, in our calculation we will not consider contributions from quark- 
disconnected diagrams. This is not a source of error given the fact that the perturbative expressions 
we will compare to will not include singlet contributions as well (entering at & (a s 3 ) for the vector 
and & (q s 2 ) for the pseudoscalar currents). But given the absence of quark-disconnected diagrams 
the two-point correlator of 7° will coincide with that of In the continuum limit vector flavor 
symmetry restoration will entail the latter to become equal to the correlation function of J^. This 
circumstance allows us to exploit the features of tmLQCD when it comes to the multiplicative 
renormalization of the bare current correlators to our advantage. 

In terms of the currents defined above the renormalized and dimensionless vector and pseu- 
doscalar moments read in the twisted basis 

G v n =Z 2 v a 6 (t/ar(4 /3 4 / \t,p = 0)r" n 

t/a=-N,/2+l 

N./2-1 

= Z 2 A a 6 £ (t/a) n (Jp^(t,p = 0)) (3.3) 

t/a=-N t /2+\ 

G p n = (2an h ) 2 ( ^) 2 a 6 (t/a) n {jf jf(t,p = 0)Y"' m 

\^PJ t/a=-N,/2+\ 
N,/2-l 

= {2a i l h fa 6 £ (t/a) n (4j^(t,p = 0)). (3.4) 

t/a=-N,/2+\ 

Using jU^ = Hh/Zp we introduced additional factors of ajih such that only the scale indepen- 
dent ratio Zs/Zp is needed for the scalar moments and no renormalization factor for the pseu- 
doscalar moments. For the scale independent renormalization factors Zp/Zs,Za,Zv we use the 
non-perturbative renormalization data provided by ETMC ([7] and private communication). 

The ensembles we choose comprise four different lattice spacings ranging from a ps 0.05 fm 
to a k, 0. 1 fm and light pseudoscalar masses in the range 280 MeV < mps < 650MeV as well as up 
to two lattice volumes. For each triple {a, L, mps) the current two-point functions were measured 
with four to seven charm valence quark masses such that the charmed meson masses anij/y/afps, 
antrjjafps and am D ±/afps in units of the light pseudoscalar decay constant covered the physical 
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value [8]. We are thus able to study the dependence of the moments G = G(a,L,mps,H c ) on all 
lattice parameters. The values of the light pseudoscalar decay constants at the physical point for all 
four lattice spacings were calculated in a separate dedicated fit along the lines of [5]. 



4. Analysis and results 
4.1 General outline 

In our analysis we will model the dependence of the moments on the lattice parameters \i q , \i c 
and a. For extrapolating to the physical light quark mass we shall use the charged pion mass mps, 
for interpolating to the physical charm quark mass the ground state mass determined from the cc 
non-singlet vector current correlator mjiy and the lattice spacing dependence will be studied using 
a/fps- Finite volume effects turn out to be negligible in the charm sector. 

We shall use two methods: 

• interpolate the lattice data at each value of the lattice spacing to common reference points 
((mps/fps) , {inj/y/fps)y ef - This strategy we shall denote with "ref". 

• perform a combined fit to our data describing the combined imps I fps, mj/y/fps, a ) depen- 
dence. This method we shall denote with "all" and it is based on splitting the fit function 
into a continuum part and one that models lattice artifacts as follows 

M N 

& (afp S , am PS , amj/ v ) = & cortt x & latt = Y,T.( amp s/ a fPs) 2 '(. am J/v/ a fPs) J 

i=0j=0 

( 

1+ £ (afp S f(amp S ) 2m (am J/v ) 2 " 

i 0<l.m,n<2 
\ 0</+m+n<2,4 



(4.1) 



We then read of the value of the moments at zero lattice spacing and at the physical point by 
setting 

G \ physical = ^cont («*//*, mj/ w /f n ) . (4.2) 

Either way we will end up with estimates for the continuum values of the moments or functions 
thereof. With these estimates we can then set up determining equations for the MS quark mass 
m c and the strong coupling a s using the perturbative representation of the moments from the low 
momentum expansion of the polarization functions. We thus set 



I physical 



£(^Yc«(m c ,A!) (4.3) 



where similar to G the coefficients C^ 1 ' will be functions of the original expansion coefficients 
in equation 2.2. In the two cases we consider we either use the charm quark mass or the strong 
coupling as input and solve the equation for the remaining quantity. Errors are estimated using a 
bootstrap method. 
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Figure 1: Example for continuum and mp$ extrapolation of G V A with the re f method. 



4.2 Moments from the vector current correlator 

The moments of the vector current correlator provide a benchmark of the method because their 
values are accessible using a dispersion integral and measurements of the hadronic cross section 
ratio R{s). We can thus make a comparison of lattice and continuum data already on the level of 
individual moments. The values of the continuum moments we compare to have been provided by 
the authors of [9] and recently in [1] (cf. the detailed description of the extraction of the charm 
piece in [9]). 

Taking into account all the explicit factors of the lattice spacing in equation 3.3 dimensional 
analysis implies the relation of dimensionless lattice moments G v n and the corresponding continuum 
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No. 


f^MKln + iy./iUn^/Ql] 1 ^ 




r Vl/(2«), 
U 2n+2 fa ^11) 


1 


0.04107 (32) 


0.04170 (25) 


0.04146 (77) 


2 


0.08792 (48) 


0.08810(52) 


0.08696 (87) 


3 


0.13081(60) 


0.13059 (68) 


0.12945(94) 


4 


0.17106 (70) 


0.17098 (82) 


0.16959(102) 



Table 1: Comparison of continuum vector moments with results from ref and all methods, 
quantities g„ at non-zero lattice spacing 

G„ = - pr _ 7 + lattice artifacts . (4.4) 

[am c ) n 1 

In figure 1 we show exemplary data for the ref method: the left panel shows the continuum ex- 
trapolation of the vector moment G\ at light pseudoscalar reference mass anips/afps = 2.5 for five 
charm meson reference masses amjj^/afps = 22.5, 23.0, . . . , 24.5 (physical point at mjiy/ ' f n = 
23.69(7) [8]). The reference points with lower light pseudoscalar masses than shown in the plot 
(anips/afps = 2.0, 2.2) are not entirely covered by the data from the coarsest lattice which is why 
we leave it out of the extrapolation and use a linear ansatz in a 2 . The right-hand side panel shows 
the extrapolation to the physical value of the light pseudoscalar mass (amps/afps) ie — > m n j ' f n = 
1.068(3) [8]. For the second extrapolation we again use a polynomial ansatz of maximally second 
degree. The "all" method gives comparable results. 

In table 1 we compare values our continuum extrapolated results for the four lowest lattice 
moments at the physical point to the continuum moments [9] (second column) determined using 
experimental data. Apart from the lowest moment M\ I G\ we find good agreement between both 
the two methods and the lattice and continuum moments. 

By comparing to perturbation theory we are now able to extract the MS charm quark mass. To 
that end we use the strong coupling as an input parameter: starting from the PDG value a s {pi = 
M z ,N f = 5) =0.1184(7) [8] we evolve it to a s (n = 3 GeV,N f = 4) = 0.255(4) using the RunDec 
program [10]. The results for the solution for the four lowest vector moments are collected in 
table 2. The first contribution to the uncertainty stems from the statistical error of the moment 
extrapolation, the physical scale {f n ) and the value of a s . The second one represents the systematic 
uncertainty from the choice of the renormalization scale: it is obtained by matching lattice and 
continuum moments at \l = (3 ± 1 ) GeV and evolving the result back to the reference scale jx = 
3GeV using 4-loop evolution [9]. 

If for each individual method we combine the quark masses from the different moments (taking 
into account their strong correlation) we find for the combined values 

o ,n f 0.979 (09) GeV (ref) 
<Mf =3GeV, W/ = 3 + l) = | 0998 ; i4 ; GeV y ■ (4.5) 

The results from both extrapolation methods turn out to be compatible with the findings of reference 
[3],m c (|i = 3GeV,iV/ = 3 + l) = 0.986(6)GeV. 

A consistency check with the lowest pseudoscalar moment using the determined charm quark 
mass as input leads to a value of the strong coupling in good agreement with the value used as input 
for the charm mass. This will be discussed in detail elsewhere. 
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No. m c (jU = 3GeV)[GeV] (ref) m c {]X = 3GeV)[GeV] (all) 
~~ 1 0.971(09) (01) 0.979 (24) (01) 

2 0.981 (10) (02) 0.998 (15) (02) 

3 0.990 (10) (11) 1.001 (12) (11) 
__4 1.014(08) (35) 1.024 (09) (34) 

Table 2: Comparison of results for the charm quark mass using re f and all extrapolated vector moments. 

5. Summary and Outlook 

With this intermediate report we showed that within the twisted mass formalism and with 
presently available statistics we can determine the moments of the charm vector current correlator 
in agreement with experimental results and with comparable uncertainty. Following two different 
analysis methods we can extract the MS charm quark mass from both methods and find agreement 
taking into account both the statistical and systematic uncertainties. Yet from the comparison of 
the central values of both analysis methods we infer that with the presently available quality of data 
a systematic error of G (20) MeV must be taken into account for the charm quark mass value. 

A consistency check with the lowest pseudoscalar moment using the determined charm quark 
mass as input leads to a value of the strong coupling in good agreement with the value used as input 
for the charm mass. Currently we also investigate other methods to extract the strong coupling from 
flavor singlet current diagrams as recently presented in ref. [11]. 

As a next step it will be very interesting to apply the methods discussed here to the Nf = 
2 + 1 + 1 gauge configurations of ETMC [12]. 

We thank all members of ETMC for the most enjoyable collaboration. This work is funded in 
part by the DFG within SFB/TR9-03. The computing time was made available to us by FZ-Julich 
on JUROPA and JUGENE. 
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